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Abstract. Using gamma-ray data from observations of the Milky Way, Andromeda (M31), and 
the cosmic background, we calculate conservative upper limits on the dark matter self-annihilation 
cross section to a number of final states, over a wide range of dark matter masses. We first constrain 
annihilation to a pair of monoenergetic gamma rays, ((7Av) rr , and show that in general our results 
are unchanged for a broader annihilation spectrum, if at least a few gamma rays are produced with 
energies within a factor of a few from the dark matter mass. We then place constraints on the 
self-annihilation cross section to an electron-positron pair, (<JAv) e + e -, usm g gamma rays produced 
via internal bremsstrahlung radiative corrections. We also place constraints on annihilation into 
the other charged leptons. We make conservative assumptions about the astrophysical inputs, and 
demonstrate how our derived bounds would be strengthened if stronger assumptions about these 
inputs are adopted. 
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INTRODUCTION 

For thermal relic Dark Matter (DM), the self-annihilation cross section must be (o^v) ~ 
3 x 10~ 26 cm 3 s _1 to obtain the observed relic abundance. There are a number of ways 
in which the self-annihilation rate can be enhanced e.g., [1, 2, 3, 4, 5, 6, 7, 8, 9], while 
non-thermal relic DM may allow larger values for (o^v) [10, 11]. Self-annihilation can 
lead to indirect detection of DM through observation of the annihilation products. We 
constrain the self-annihilation cross section by considering the gamma rays produced 
via DM annihilation to a number of final states. 

We first focus on annihilation to a pair of monenergetic gamma rays. This annihilation 
channel has a very clean signature, with Ey = m x , and future observation of such a 
line could provide clear evidence for the DM mass. We compare the calculated signal 
with gamma-ray observations of the Galactic Center, Andromeda (M31), and the cosmic 
background, and place upper limits on the partial self-annihilation cross section, (oxv)yy- 
Our limits cover a broad DM mass range, from 20 keV up to 10 TeV. We show our 
constraints are applicable even if DM annihilates to a broader spectrum of gamma rays, 
in general constraining the total annihilation cross section. With a conservative choice 
for the branching ratio, Br(yy) = 10~ 4 , we constrain the total annihilation cross section, 
(OAv)totai = (o A v) rr /Br('YY). 

Next we place constraints on the self-annihilation cross section to an electron-positron 
pair, {<JAv) e +e- ■ These particles will inevitably be accompanied by gamma rays through 
a process known as internal bremsstrahlung (IB) ox final state radiation. This is the 
emission of a gamma-ray from one of the final state charged particles, occurring at the 



Feynman diagram level. See Refs. [12, 13, 14, 15] for a detailed discussion. We use 
these gamma rays to constrain (<JAv) e + e - by again comparing the calculated annihila- 
tion signal with the observed gamma-ray flux. We place corresponding constraints on 
annihilation to the other charged leptons. The IB gamma-ray spectrum per annihilation 
does not depend on the tree level cross section, so our analysis remains DM model in- 
dependent. Because emission occurs at the Feynman diagram level, the gamma-ray flux 
from IB remains independent of uncertain astrophysical phenomenon such as the inter- 
stellar magnetic field or radiation field. 

The e + e~ final state is especially interesting in light of the recent PAMELA/ ATIC/ 
Fermi results suggesting an excess in the observed positron flux. There are now nu- 
merous DM models that can explain these results, many with a large branching ratio to 
charged leptons, e.g. [3, 4, 5, 6, 7, 8, 9]. Our limits on (oAv) e + e - directly constrain many 
models attempting to explain the positron excess. 

We show how our results depend on the astrophysical uncertainties in the DM den- 
sity profile. We are conservative in our analysis methods and input choices, and show 
how our results would be strengthened by more optimistic choices. We compare our 
constraints with pre-existing bounds on the total annihilation cross section. We find that 
our constraints on (oxvKotai are particularly strong for small DM masses, while our con- 
straints on {gav)yy are stronger than those on (oAv) e + e - by a factor of < 10 2 , as expected 
from the a ~ 1/137 suppression of the radiative corrections. 



ANNIHILATION TO yy 

Our analysis and results are described in detail in Mack et al. [16]. The flux of gamma 
rays from DM annihilation in our halo is given by 

d<Z> y _ (o A v) 1 dNy 

dE 2 Jo Anm 2 % dE ' 

where dN y /dE is the gamma-ray spectrum per annihilation, AO. is the field of view in 
steradians, and 

2% fV f e ~ max o / / \ 

js£i = J ° y y p ( v ^ c ~ 2iR&c cos ¥ + i2 ) di sin wd¥ ' (2) 

is proportional to the dark matter density squared, integrated over the line of sight. For 
annihilation to a gamma-ray pair, the spectrum per annihilation is simply dNy/dE = 
28(m x -E). 

We collect gamma-ray flux data from a variety of sources, covering a large range of 
energies. We compare the calculated annihilation flux with the observed flux, integrating 
each over an energy bin of equal size, and find an upper limit on (<7av) by demanding 
that the annihilation flux be smaller than the entire observed flux. This is extremely 
conservative, as in reality, the DM annihilation flux would only be a small fraction of 
the observed flux. 

We focus on the Galactic Center flux, using data from COMPTEL and EGRET aboard 
the CGRO [22, 23], INTEGRAL [21], and H.E.S.S. [24]. We also look at gamma rays 
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FIGURE 1. The limits on the partial cross section, (oAv)yy, derived from the various gamma-ray data. 
Our overall limit is shown as the dark shaded exclusion region. For comparison, the light-shaded region 
shows the corresponding limits for the NFW (rather than the Kravtsov) profile. 



arising from annihilation in the M3 1 (Andromeda) galaxy, and from cosmic annihilation, 
using the methods described in Mack et al. [16]. We use M31 data from EGRET, 
CELESTE [25] and HEGRA [26], and diffuse cosmic flux data from INTEGRAL, 
COMPTEL, EGRET and the Solar Maximum Mission [27]. 

We place constraints on (oav)yy = (°Av)total-^ r (T7)> f° r both the conservative 
Kravtsov [18] and the commonly-adopted NFW DM density profile [19], over an 
energy range spanning 20 keV to 10 TeV. These results are shown in Fig. 1. 

To find a constraint on the total cross section, we need the branching ratio to this 
final state. Typical branching ratios are ~ 10~ 3 . In Fig. 2 we show our upper bound 
on (cTAv)totai using a conservative choice of Br(yy) = 10~ 4 . In this plot we use the 
conservative Kravtsov density profile. For comparison, we show other limits on (o^v), 
including a bound based on observations of the Neutrino flux, from Yiiksel et al. [17], 
and the standard cross section for thermal relic DM. The KKT and Unitarity bounds are 
described in Ref. [16]. 

Most of these experiments had modest energy resolution. To be conservative, we 
assume an analysis bin with a logarithmic energy width of 0.4 in \og l0 E (i.e., A(ln£') ~ 
1) for the Galactic and cosmic diffuse analyses; this is at least as wide as the energy 
bins reported by the experiments. In effect, our results are what one would obtain for 
an annihilation gamma-ray spectrum as wide as 0.4 in \og lQ E. The exception is the 
INTEGRAL line search. 
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FIGURE 2. The gamma-ray and neutrino limits on the total annihilation cross section, selecting 
Br(yy) = 10~ 4 as a conservative value. The unitarity and KKT bounds are also shown. The overall bound 
on the total cross section at a given mass is determined by the strongest of the various upper limits. 



ANNIHILATION TO e + e~ 



For annihilation to an e + e pair, we follow a similar method to that outlined above. 
We constrain (<JAv) e + e - by demanding that the calculated gamma-ray signal from IB 
emission not exceed the entire observed gamma-ray flux. Our analysis and results 
are described in detail in Bell and Jacques [12]. The IB gamma-ray spectrum per 
annihilation is 



dN v 



1 dom 



dE Otot dEy 



where 



dom oc 
dE 



r 



EyTt 



In 



4m x (m x 



m-, 



1 + 



4m x (m x 



4m| 



(3) 



(4) 



and Otot lS the tree-level cross section for annihilation to e + e~ . Note that o to t factors 
out from the IB cross-section. This important feature implies that the IB spectrum 
is independent of the unknown physics which mediates the lowest order annihilation 
process. Fig. 3 shows E 2 dN r /dE for IB emission for a number of DM masses. 

We use Galactic gamma-ray data from COMPTEL, EGRET and H.E.S.S. We also use 
data from CELESTE observations of M3 1 to cover the gap between the highest energy 
EGRET data and the lowest H.E.S.S. data points. The strongest bounds are obtained 
by choosing an analysis bin where the ratio of the calculated annihilation signal to the 
observed background is largest. Fig. 3 shows that dNy/dE falls with energy. However, 
the observed gamma-ray flux falls with energy at a steeper rate, so that the strongest 
results are obtained by choosing an analysis bin with the DM mass as the upper energy 
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FIGURE 3. Internal bremsstrahlung gamma-ray spectra per XX ~ * e + e~ annihilation, for m z = 100 
GeV, 200 GeV, 500 GeV, 1000 GeV. 
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FIGURE 4. Upper limit on (ov) e + e - as a function of DM mass for the Kravtsov (solid), NFW (dashed) 
and Moore (dot-dashed) profiles. 



limit. In Fig. 4 we show the upper limits on (<JAv) e + e - as a function of DM mass, 
using the observational data described above. We give the Galactic Center results for 
the Kravtsov and NFW profiles mentioned earlier, as well as the steep Moore profile 
[20]. We show the CELESTE constraint using only the Kravtsov profile. 
In Fig. 5 we show the upper bounds on the annihilation cross sections into e + e ~, 
, and t + t~, based upon the IB emission from each final state (all use the conser- 
vative Kravtsov profile). For annihilation to t + t~, we do not take into account gamma 
rays arising from the hadronic decay modes. Including this would strengthen our bounds. 
For comparison, we also show limits on {gav)yy discussed in the previous section, as 
well as some earlier bounds on (o^v). For full details see Ref. [12]. 



10 -3 -2 -1 1 2 3 4 5 

10 10 10 10 10 10 10 10 10 
m x [GeV] 

FIGURE 5. Upper limits on Br(ii) x (<Jv) totu i for various final states ii — e + e~ (solid black; labelled), 
(thick dashed; labelled), T + T~ (thick dashed; labelled), yy (solid; labelled), and vv (solid; la- 
belled), using the conservative Kravtsov profile. Also shown are the KKT (thin dashed) and unitarity (thin 
dot-dashed) limits on the total cross section described in Refs. [16, 12], and the cross section for thermal 
relic DM (natural scale). The vv limits are taken from Ref. [17]. 



DISCUSSION 

While the Kravtsov, NFW and Moore profiles diverge towards the center of the Galaxy, 
they are similar at large radii. As the INTEGRAL, EGRET and COMPTEL observations 
encompass relatively large angular scales, the density profile changes have a modest 
effect. On the other hand, the H.E.S.S. constraints correspond to a much smaller angular 
region toward the Galactic Center, and vary by orders of magnitude depending on the 
profile adopted. (See Ref. [17] for a full discussion of the differences between the 
profiles for different angular regions.) 

If DM annihilates to e ± , photons will be produced not only by IB, but also by 
energy loss processes including inverse Compton scattering and synchrotron radiation. 
However, the intensity of these signals depend on uncertain astrophysical parameters, 
such as magnetic field strength, radiation field intensities, and electron diffusion scales. 
By contrast, IB is free of these astrophysical uncertainties, and has a fixed spectrum 
and normalization. Another key difference is the energy of the photons. Synchrotron 
radiation produces generally low energy photons, while IB provides some hard gamma 
rays near the endpoint. Since the background flux falls off with energy, these hard 
gamma rays are extremely useful. The sharp edge in the IB spectrum at E = m x can 
be used to diagnose the DM mass; this is not possible with synchrotron radiation. 

Refs. [3, 4, 5, 6, 7, 8, 9] and others have recently proposed models in which DM 
annihilates directly to charged leptons, with cross sections well above that expected for 



a thermal relic. This may account for anomalies in cosmic ray spectra from PAMELA, 
ATIC and Fermi, and microwave signals from WMAP, all of which seem to require more 
electrons and positrons than can be explained otherwise. Our bounds on {cav)i+i- will 
directly constrain the allowed parameter space for these types of DM models. 
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